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The Blume-Capel model with infinite-range interactions presents analytical solutions in both 
canonical and microcanonical ensembles and therefore, its phase diagram is known in both ensembles. 
This model exhibits nonequivalent solutions and the microcanonical thermodynamical features 
present peculiar behaviors like nonconcave entropy, negative specific heat, and a jump in the 
thermodynamical temperature. Examples of nonequivalent ensembles are in general related to 
systems with long-range interactions that undergo canonical first-order phase transitions. Recently, 
the extended gaussian ensemble (EGE) solution was obtained for this model. The gaussian ensemble 
, and its extended version can be considered as a regularization of the microcanonical ensemble. 

■ They are known to play the role of an interpolating ensemble between the microcanonical and the 

canonical ones. Here, we explicitly show how the microcanonical energy equilibrium states related 
to the metastable and unstable canonical solutions for the Blume-Capel model are recovered from 
EGE, which presents a concave "extended" entropy as a function of energy. 
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4-> ■ I. INTRODUCTION 

> 

' Thermodynamical properties of physical systems are described in Statistical Mechanics by the ensemble theory, 
"j^ Its development rests on the existence of both additivity and thermodynamic limit, which are essential ingredients 
to describe critical phenomena. The existence of a thermodynamic limit is related to the extensivity property which 
is required in order to obtain meaningful physical results in that limit Large enough systems with short-range 
'■^ ■ interactions present the additivity property for extensive quantities like energy and entropy. Additivity means that 
^ ^ when we consider the energy, or any other extensive quantity of a system, and divide the system into two or more 
O ,: macroscopic parts, the total energy is still the sum of the energies of those parts. From theoretical backgrounds, 
the canonical and grand-canonical ensembles approximate the microcanonical ensemble in the limit of infinitely large 
number of particles, where surface effects and fluctuations can be disregarded with respect to the bulk mean values 
0,0]. There is definitely more consensus around the microcanonical ensemble as the most fundamental one. Thus, 
[ the canonical and grand-canonical ensembles are taken as approximations for the microcanonical ensemble whenever 
the full equivalence does not hold. In cases where full consistency of statistical ensembles holds for systems that 
undergo phase transitions, it is found that finite-size scaling relations still place the microcanonical approach as the 
fundamental one Q]. 

There are many examples of systems whose microcanonical equilibrium properties contrast to the ones obtained 
in the canonical ensemble. Differences in the equilibrium properties have been described analytically for systems 
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' with long -range interactions [5|-|8j , which are consequences of the nonconcavity of the entropy as a function of energy 
0, B lid ] The nonequivalence appears for those systems when the canonical ensemble exhibits a first-order phase 
transition and this has been demonstrated by means of the large deviation theory [llj . Such nonequivalence results in 
uncommon features like temperature discontinuity and negative specific heat in the microcanonical ensemble [l0l-[l3| . 
In this paper, we will be dealing with a less popular ensemble, the extended gaussian ensemble (EGE) [1^|, which 
5—1 ' is constructed upon the gaussian ensemble |16l42C| . They are known to present a smooth interpolation between their 
limiting behaviors, corresponding to the microcanonical and canonical ensembles. Recently, the EGE solution was 
obtained for the Blume-Capel (BC) model with infinite-range interactions [2ll|. Our numerical calculations illustrate 
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this limiting behavior in a model with a tricritical point. This tricritical point does not coincide numerically in the 
two standard ensembles. Here, we explicitly show how the microcanonical energy equilibrium states related to the 
metastable and unstable canonical solutions for the Blume-Capel model are recovered from EGE, which presents 
a concave "extended" entropy. Moreover, the discontinuous microcanonical solution is also investigated and the 
metastable and unstable lines are also obtained from EGE. 

The content of our presentation is as follows. The EGE analytical solution for the BC model is reviewed in the next 
section. Analytical and numerical results for thermodynamical quantities characterizing this ensemble are presented 
in section 3. In section 4 we summarize the results of this methodology in identifying the region of nonequivalence of 
ensembles only from the EGE point of view. 



II. EXTENDED GAUSSIAN ENSEMBLE AND THE BLUME-CAPEL MODEL 



The EGE is defined by the condition 7^0 and probability density 

where Zj{U, a) stands for the normalization constant, which is the corresponding partition function in EGE (l5l[22l . [23l |. 
As usual, p{E) is the density of states. Here, 7 and a are parameters, and U stands for the mean energy, which is a 
and 7 dependent, U = U{a,^). The extended gaussian ensemble is a particular case in a class of ensembles defined 
by functions g{E) [l^, H^, where the quadratic form g{E) = j{E — UY is a convenient choice. The mean energy U 
can be obtained from the extended Gibbs's measure, 

U = j EP^,a{E)dE. (2) 

This ensemble has two important limits. The canonical results correspond to 7 ^ and in this case the parameter 
a identifies the inverse thermodynamical temperature, a = l/ksT. On the other hand, the limit 7 ^ 00 corresponds 
to the microcanonical case. This can be seen as \\m^f^ao\/T^ H Z~f{U ^ a) through the use of the Dirac's delta sequence 
in the gaussian form (23 | . The finite 7 case presents an intermediary thermal description between the known limiting 
ensembles. All thermodynamical quantities can be defined in this ensemble as well. The extended thermodynamical 
potential is analogously defined to the canonical ensemble, (£/,«) = —\x\Z^{U^a). From here, the derivative at 
fixed value 7 gives the mean energy, 
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which parallels that of the usual canonical approach. The extended heat capacity is [iHli^ ^7 = —a^{dU /da)^. 

Let us now consider the Blume-Capel model [2^ [2^, a particular case of the Blume-Emery-Griffiths model [27 
Here, we work with its infinite-range interaction version, 

N 



4=1 




H{S)^Aj:sf--L^ Y.^^ ' (4) 



where Si = {0,±1}. The couplings J > and A are the exchange and crystal-field interactions, respectively. This 
is a well studied model, which presents a rich phase diagram in the (A/J, T/J) plane. It exhibits a first-order phase 
transition line, tricritical point, and a second-order phase transition line. The critical pro per ties of the BC model 
have been studied analytically in both the microcanonical Q and canonical ensembles pvL l28j|. The BC model has 
become an outstanding model for what concerns the nonequivalence of ensembles. The canonical and microcanonical 
ensembles do not yield the same phase diagram for the first-order critical line 5\. The canonical tricritical point 
occurs at (A/J, T/ J) ~ (0.46209812, 1/3), which gives origin to the canonical first-order phase transition line for larger 
values of A/J (and lower values of T/J). The microcanonical solution identifies the tricritical point at (A/J, T/J) ~ 
(0.46240788,0.33034383). 

The EGE analytical solution [2l| is obtained from the extended partition function by means of Hubbard-Stratonovich 
transformations. To this end, we introduce the order parameters magnetization M = X]i=i ~ ^+ ~ the 
quadrupole moment Q = X^^Li = + iV_, where Af+ and are, respectively, the number of sites with up and 
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down spins. Also, let Nq be the total number of zero spins, then N = 7V+ + + Nq is the total number of spins in 
the system. Therefore, the extended partition function reads 

N N-No . 
No=0 N+=0 +' 

Here, it is convenient to work with the intensive quantities q — Q/N and m = M/N . Let us also define K = J/2 A 
and £ = U/AN the energy per site, following the authors in [5]. 

First of all, we recall the EGE solution for 7 — > cx), which yields the known microcanonical ensemble results. The 
extended partition function must be well-behaved in this limit. Thus, it is required in equation ([5]) that e = q — Km? , 
which becomes a constraint equation for the mean energy e as 7 — > 00. The thermodynamic limit N ^ 00 can be 
evaluated and one obtains the thermodynamical potential per site 



<p(e, a, m) ~ eaA 



(6) 



as a function of mean energy per site e and, parameters a and mean magnetization m. The term inside the brackets 
is the correct microcanonical entropy — Smicro(£, w) obtained in Q as a function of e and to. The entropy is an 
even function of m and a nonconcave function of independent variables to and e. The nonconcave behavior is the 
necessary condition for metastability and instability in the canonical states [2^ and gives origin to the nonequivalence 
of ensembles (microcanonical and canonical). The consequences on the thermodynamical quantities specific heat and 
specific susceptibility have been analyzed in [2l| . 

The solution for finite 7 describes up to what extension one recovers the microcanonical stable states. The full 
equivalence of these results with the stable states of the microcanonical ensemble is achieved for finite 7 only for A/ J 
between the canonical and the microcanonical tricritical points. For couplings A/J larger than the microcanonical 
tricritical point, i.e., on the microcanonical first-order transition line, one needs 7 — > (X) for such full recovery of the 
microcanonical states [2ll |. 

The finite 7 solution for the thermodynamical potential per site reads 



TO , / a + TO 



^,(£,a,TO,g) = ,lnl^^^— ^l+:^ln(^^j+ln(l-g) 



The EGE analytical solution poses a clear way of studying the formation of stable energy-dependent microcanonical 
states as the parameter 7 increases. Furthermore, it is possible to draw the EGE solutions to observe the discontinuous 
microcanonical transition. This will be illustrated in the next section. 



III. RESULTS AND DISCUSSION 

The lack of ensemble equivalence is a consequence of a nonconcave entropy. As discussed in [2^, the possible 
stable states, labeled by e, are the ones defined by the microcanonical entropy Sinicio(e)- If those values of s are 
not global minima of the canonical thermodynamical potential per site as we vary the temperature, then we do not 
find such correspondence between these states in both ensembles. Local minima and local maxima in the canonical 
thermodynamical potential per site are, respectively, metastable and unstable states. 

To see how the nonconcave entropy for the BC model induces this nonequivalence of ensembles, we initially evaluate 
the caloric curve T versus e, in the canonical limit. To this end, we calculate the equilibrium solutions of ifijis, a, to, q) 
for 7 — >■ 0. The equilibrium solutions are easily obtained in this limit and are shown in figures 1 and 2. Figure 1 shows 
caloric curves for different values of A/J. Notice that the point A/J = 0.462098, although close to the canonical 
tricritical point, is still in the second-order phase transition region according to both canonical and microcanonical 
predictions. The points A/J = 0.4622 and 0.4623, are in the canonical first-order phase transition region, but they 
are in the second-order region from the microcanonical point of view. Figure 2 shows caloric curves for A/J in a 
region where both ensembles present first-order phase transition lines. Notice that we have not drawn in those figures 
horizontal lines corresponding to the canonical Maxwell construction. Such constructions can be appreciated in (5| for 
some values of A/ J. The missed stable microcanonical states are the ones related to the gaps in e in both figures. We 
anticipate that the presentation of those stable microcanonical states is left to figure 4, while in figure 5 we included 
all the equilibrium solutions. 
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FIG. 1: EGE temperatures obtained in the canonical 
limit for some values of A/ J. A/J = 0.462098 is in 
both canonical and microcanonical second-order phase 
transition regions. A/J = 0.4622 and 0.4623 are in the 
canonical first-order phase transition region but in the 
microcanonical second-order ( — 
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FIG. 2: EGE temperatures obtained in the canonical 
limit for some values of A/J. The values of A/J are 
in both canonical and microcanonical first-order phase 
transition regions. 




The EGE equilibrium states for finite 7 can be obtained by solving the extremum equations for Lp^{e,a,m,q). 
We call attention to the fact that the microcanonical constraint for the specific quantities is not enforced here, the 
variables £,to and q are treated as independent variables in this approach. The linear term in 7, equation ([7]), can be 
seen as a constrained equation leading to the microcanonical ensemble only for 7 — >■ cx). 

As done in (2l| . the desired solutions must satisfy equation ([3]), dip-^/da = eA, and the extremum conditions 
dif-f/dm = 0, and d(£^/dq — 0, for a fixed e. The stability criterion in the canonical ensemble for the BC model is 
nicely discussed in We must observe that the solutions for the three equations above are not stable for all 7. 
Since the analytical expression ipj comes from a saddle-point approximation, one needs to study the stability of those 
EGE solutions as a function of m and q. To this end, the determinant of the Hessian matrix, 

S IS ' 

dqdm dq^ / 

is analyzed in the T versus e plane as a function of 7. This amounts to exploring which points {m, g} minimize ip-^ 
for fixed T and e, and satisfy the condition dlrn, q) > 0. 

Figure 3 collects our results of stable points for different values of 7. That behavior is illustrated for A/ J = 0.462407, 
a coupling close to the microcanonical tricritical point and on the canonical first-order critical line, but still on the 
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FIG. 4: EGE temperatures obtained in the micro- 
canonical limit for A/J = 0.4622 and 0.4623. The 
couplings A/J are in the microcanonical second-order 
phase transition region. 
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FIG. 5: EGE temperatures obtained in the micro- 
canonical limit for A/J = 0.4625 and 0.4627. The 
couplings A/J are in the microcanonical first-order 
phase transition region. 



second-order critical line according to the microcanonical calculations. This figure shows that, as 7 increases, one 
recovers the microcanonical solution. It can be shown that for sufficiently large 7, s^{e) becomes entirely concave and 
continuous on e [2I, □01 , 



2^ 



5e2 



< 0. (8) 



In other words, the addition of the term in 7 to the thermodynamical potential changes the energy range where its 
nonconcavity is observed. This energy range is reduced as 7 increases. Equation ([5]) can be used to find the necessary 
condition on 7 to obtain the microcanonical solution, 

^ ^ 2Ar2(e)c(£)' 

Here, c(e) is the specific heat, a negative quantity for energy ranges in figure 3 such that dT/de < 0. Figure 3 shows 
that, for the coupling A/J — 0.462407, a finite value of 7 is able to recover all stable microcanonical states. Actually, 
7 ~ 4950 is the minimum value necessary for this end. To obtain the corresponding stable microcanonical states 
for A/J = 0.4622 and 0.4623 in figure 1, one needs a rather small value for 7. It can be shown 21] that 7 30 is 
sufficient for such couplings A/J. The numerical solutions are presented in figure 4. 

For A/J larger than 0.46240788, the microcanonical tricritical point, the recovery of all microcanonical states is 
achieved only for 7 — > cx). Different of a continuous recovery of microcanonical states, as shown in figures 3 and 4, the 
caloric curves on this microcanonical first-order phase transition line sets a new physical situation. Here, one observes 
jumps in the temperature, the analogous of the latent heat in the canonical ensemble. The energy where this transition 
occurs can be obtained by defining a vertical line corresponding to the Maxwell construction for the microcanonical 
ensemble [l^[3l|. Moreover, the microcanonical hysteresis effect is also present. This effect in the thermodynamics of 



a self-gravitating system has been discussed by Chavanis [3lJ. Our figure 5 illustrates this situation for A/J = 0.4625 



and 0.4627. Although for such values of A/J the full recovery of all microcanonical states is only possible in the 
limit 7 ^ 00, we can still obtain valuable estimates about those states by choosing a very large numerical value for 
7. Moreover, we were able to draw the continuation lines that define two branches, the high and low temperature 
branches, which produce the unstable and metastable lines of microcanonical states. The solutions depicted in figure 
5 show that T(s) is a multivalued function of e, signaling a microcanonical first-order phase transition, in addition 
to the standard profile of the canonical first-order phase transition when one draws horizontal lines corresponding to 
the Maxwell's construction. 



IV. CONCLUSIONS 



The EGE approach leads to analytical expressions for the extended free energy and entropy in a simple way, and 
quantifies the nonequivalence of ensembles. The nonconcavity of the microcanonical entropy function means that the 
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system contains several energy-dependent stable states, revealed in the microcanonical ensemble, which have their 
counterpart in the canonical ensemble, but they are unstable or metastable in this latter ensemble. By microcanonical 
entropy we do not refer to the standard Legendre transform of the canonical thermodynamical potential, which always 
produces a concave function of energy. The new term in 7 added to the canonical free energy turns such unstable 
or metastable points into equilibrium points in the extended ensemble i22i . i23i l29l [s^ because it presents a concave 
"extended" entropy. Moreover, it is explicitly shown for the BC model that not only the stable states present in the 
microcanonical approach can be obtained from EGE but also the unstable and metastable ones. This is illustrated in 
figure 5, which shows the complexity of the microcanonical phase diagram and the usefulness of EGE for such study. 
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